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Abstract 

Little Higgs models with T parity contain an attractive dark matter candidate, the 
heavy photon. We compute the cross section of the heavy photon annihilation into Z- 
photon pairs, which turns out to be substantially higher than the previously computed 
cross section for the two photon final state. Unfortunately, even with this enhancement, 
the monochromatic photon flux from galactic heavy photon annihilation is unlikely to 
be detectable by GLAST or the currently operating atmospheric Cerenkov telescopes. 
We also compute the flux of high-energy neutrinos from the annihilation of the heavy 
photons captured by the Sun and the Earth. The maximum flux of upward-going 
muons due to such neutrinos is about 1 yr _1 km -2 . 



1 Introduction 



The presence of dark matter in the universe has been firmly established by observations of 
galaxy rotation curves, large scale structure, and cosmic microwave background radiation. 
The microscopic nature of dark matter, however, remains unknown. One attractive sce- 
nario is that it consists of stable weakly- interacting massive particles (WIMPs), with masses 
around the weak scale, ~100-1000 GeV. If the WIMPs were in thermal equilibrium with 
other species in the early universe, their relic abundance naturally matches the observed 
dark matter density. Further, many theoretical extensions of the Standard Model (SM) 
of particle physics include such particles. The most famous example is no doubt the neu- 
tralino of most SuperSymmetric (SUSY) models; another well known example is the lightest 
Kaluza-Klein excitation in models with universal extra dimensions. 

In this paper, we consider the dark matter candidate emerging in another popular class 
of theories, the Little Higgs models [1] (for reviews, see [2, 3]). In these theories, the Higgs is 
a pseudo-Nambu-Goldstone boson associated with spontaneous global symmetry breaking in 
an extended electroweak sector, which occurs at a scale / ~ 1 TeV. The Higgs acquires a mass 
due to an explicit breaking of the global symmetries by gauge and Yukawa interactions, but 
a special "collective" manner in which this explicit breaking occurs ensures that the Higgs 
is relatively light, m h <C f. Several models based on this idea have been constructed, most 
of which are tightly constrained by precision electroweak data. However, introducing an 
additional discrete symmetry, the T parity, weakens the precision electroweak constraints 
significantly [4]. The Littlest Higgs model with T parity (LHT) [5] provides an example 
of a fully realistic and natural model consistent with precision electroweak data [6], and is 
sufficiently simple to allow for detailed phenomenological analyses [7, 8, 9, 10]. 

An interesting consequence of T parity is the stability of the lightest T-odd particle 
(LTP), typically the heavy partner of the hypercharge gauge boson B H . (This particle is 
often referred to as the "heavy photon", even though it does not couple to the electric 
charge.) The Bh is weakly interacting, has a mass in the range 80 — 500 GeV, and it was 
shown to have the correct relic abundance to account for all of the observed dark matter 
in certain regions of the parameter space [7, 11]. Potential non-gravitational signatures of 
the heavy photon dark matter have been discussed in [11] (direct detection and indirect 
detection via high-energy photons produced by Bh annihilations in the galactic center) 
and in [12] (indirect detection via high-energy positrons). In this paper, we continue the 
evaluation of indirect signatures for this dark matter candidate. First, we complete the 
previous analysis of the monochromatic photon flux [11] by computing the flux from the 
B H B H — > Z^i reaction, whose cross section turns out to be substantially larger than the 
BhBh — > 27 reaction analyzed in [11]. Unfortunately, even with this enhancement the 
near-future prospects for the detection of the gamma line do not look promising, especially 
when the EGRET constraint on the continuous photon flux is taken into account. Second, 
we evaluate the neutrino flux from the annihilation of the heavy photons captured by the 
Sun and the Earth. While the predicted fluxes are too low to be observed by the near-future 
detectors such as ICECUBE, a further improvement of two-three orders of magnitude in 
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detector sensitivity would allow to probe interesting parts of the parameter space. 

2 The Model 

The LHT model has been described in detail elsewhere [5, 6, 7]; here, we briefly summarize 
the features of the model important for our analysis. The extended electroweak sector of 
the LHT model has a global 577(5) symmetry, which is spontaneously broken to 5*0(5) at 
a scale / ~ 1 TeV. At energies below the cutoff A ~ Anf, the dynamics of this sector is 
described by a non-linear sigma model. To incorporate gauge interactions, the subgroup 
Q = (SU(2) x U(l)) 2 of the SU(5) is weakly gauged. The T partity interchanges the two 
577(2) x U(l) factors. At the scale /, the gauge symmetry Q is broken down to SU (2) x £7(1), 
identified with the SM electroweak group. The gauge bosons corresponding to the broken 
generators, Wjj, W H and B H) are T-odd and acquire masses at the scale f: 1 

M(W H )~gf, M = M(B H ) ^-J= « 0.16/. (1) 

The uneaten SU(5)/ 50(5) Nambu-Goldstone bosons decompose into a T-even SU(2) dou- 
blet H, identified with the SM Higgs, and a T-odd SU (2) triplet $, which acquires a mass 
= \f2rrihf jv at one loop. In the fermion sector, the LHT model contains vector-like 
T-odd partners for the left-handed SM quarks and leptons. We will assume a common mass 
scale M for all of these particles. This choice is motivated both by simplicity and (in the case 
of the T-odd quarks) by flavor constraints [13]. Finally, the model contains an additional 
pair of weak-singlet, charge-2/3 quarks: T + (T-even) and T_ (T-odd). These are required 
to ensure the cancellation of the one-loop quadratic divergence in the Higgs mass parameter 
from the SM top. The mass M T _ is a free parameter, and 

/ 2f2 \-i/2 

With our assumptions, the spectrum and the (renormalizable) couplings of the LHT model 
are completely described in terms of the three parameters, /, M, and M T _, in addition to 
the familiar parameters of the SM. 

Due to the smallness of g' and the favorable group theory factor, the "heavy photon" 
Bh is substantially lighter than /, and is always the lightest among the T-odd bosons of 
the theory. The fermion mass parameters M and M T _ are at the scale /, so it is natural to 
consider the part of the parameter space where the Bh is the lighest T-odd particle (LTP). 
In this region, the heavy photon is stable and can contribute to dark matter. In this paper 
we assume that all of the observed dark matter density is due to the Bh. Calculations of 
the relic density have shown [7, 11] that there are three regions of the parameter space in 
which this is possible. In the first two, the B H density in the early universe is controlled by 

1 In Eq. (1) and throughout this paper, we neglect corrections of order v 2 /f 2 . 
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Bh pair-annihilation through an s-channel Higgs. The correct relic density is obtained when 
the center of mass energy of a Bh Bh collision in the nonrelativistic regime is near, but 
slightly displaced from, the resonance. There are typically two solutions, one on either side 
of the resonance. In the third region, the Bh relic density is set by coannihilation processes 
with T-odd fermions, and hence M must be close to M. The precise measurement of the 
present dark matter density implies tight correlations between the model parameters in these 
regions. It was found in [11] that they can be approximately described by 



M = < 



(m^/2.38) — 10 "Low" pair annihilation region 

(mh/ 1.89) + 44 "High" pair annihilation region (3) 

M — 20 coannihilation region 



(All masses are in GeV.) With the reduced uncertainties from the 3-year WMAP data set [14], 
n dm /i 2 = 0.104 ± 0.009, the allowed variation of M around the central values in (3) is at 
most of order ±5 GeV at the 2a level. We will thus treat the allowed regions as lines defined 
by Eq. (3). 



3 Monochromatic Photon Flux from BhBh~^ Zj/hj 
Annihilation 

A promising avenue for indirect WIMP detection is through high energy gamma rays [15]. 
Since the galaxy is transparent to photons in the interesting energy range, any features in 
the photon spectrum are preserved. This can both aid in distinguishing the signal from 
the background, and provide information about the WIMP properties. In particular, a pair- 
annihilation of WIMPs into two-body final states containing photons leads to a flux of nearly 
monochromatic gamma rays, allowing for efficient background subtraction. In the context 
of the LHT model, the monochromatic gamma ray flux from the process BhBh^ 27 was 
analyzed in Ref. [11]. In this section, we will analyze the additional two processes which 
produce monochromatic photons, BhBh^ or hj. 

The BhBh annihiliation in the LHT model is dominated by processes with an s-channel 
Higgs boson exchange. The cross section for producing a given final state X in this channel 
can be related to the partial decay width of an (off-shell) Higgs [11]: 

<A 2 s 2 - 4sM 2 + 12M 4 T(h^X) 

a x u = a(B H B H ^ X)u = ——- -5- — -= . (4) 

72M (s - m 2 h ) + m 2 h T 2 h Vs 

Here, u is the relative velocity of the annihilating WIMPs; s m 4M 2 in the galactic centre; 
and the hat on T indicates that the substitution m h — > y/s should be made in the expression 
for the on-shell Higgs decay width. 

Using the well-known results for the partial Higgs decay width in the Z7 channel [16, 17], 
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we obtain 2 



a g 2 s 3 / 2 
2048tt 4 m 2 w 



A F + A G 



1 - 



mi 



(5) 



Here, ^4jc and Ag are the contributions to the matrix element from fermions and gauge 
bosons, respectively: 



A F 
Ag 



-4^^L^N Cf Q f V f [h (r f , A/) - I 2 (r f , A/)] , 
/ 9 T f 



Et [ %( 4 (3-»^) + 



1 + - f 



5 + ^ 



'<?/ J 



A(r 9 ,A 9 ) ,(6) 



where — 0.548 is the tangent of the Weinberg angle, and 

Am 2 . Am 2 



(7) 



The sums run over all particles of the relevant type (including both the SM and the additional 
particles of the LHT model), with the electric charge, the multiplicity and the coupling to the 
Z boson of each particle given by Qi, N Ci , and V$, respectively. For fermions, V* denotes the 
vector part of the HZ coupling; there is no contribution from the axial part. In particular, 
for the extra vector-like fermions of the LHT model, we obtain 



Vi = -?- (r s (<) - s 2 w Q(i)) , 

Cw v y 



(8) 



where T 3 (U,N) = 1/2, T 3 (D,E) = -1/2, and T 3 (T+, T_) = 0. The fermion and vector boson 
trilinear couplings to the Higgs are given by yf/y/2 and CigMy/ff' 1 ' ', respectively. (With this 
normalization, i/f are the standard Yukawa couplings for the SM fermions and = 1 for 
the SM W ± boson.) The functions J lj2 are given by 



AM) 

I 2 (a,b) 



a b 



2u2 



+ 



a 2 b 



2(a-b) 2(a-b) 2 
-2(^) [/(a) - /(6)] ' 



If (a) - f(b)} + 



a 2 b 



[a - by 



[9(a) -g(b)} , 



where 



m 

g(x) = 



sm 

_i 

4 



log - 



VK 



\ix > 1 
if x < 1 

ifx > 1 
ifx < 1 



(9) 

(10) 
(11) 



2 We neglect the contribution to the Higgs width from the loops of the T-odd scalar $. 
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Figure 1: Photon fluxes from WIMP annihilation into Z^ (solid line) and 77 (dashed line), 
in the "high" pair-annihilation region of the LHT model. The fluxes in the "low" pair- 
annihilation and coannihilation regions are similar. The fluxes are normalized to JAQ = 1, 
and scale linearly with this parameter. 



For a telescope with line of sight parameterized by ^ = {6, 0) and an angular acceptance 
Af2, the anomalous photon flux due to B H B H ^ is given by [15] 

• = (5.5X 10- cm-) (0} (12^X) 2 J(*,A fi )Af ! , (12) 

where the function J contains all information about the dark matter distribution in the halo. 
The photon energy is 

E - - M t 1 " • < 13 > 

The thermal broadening of the line is much smaller than the energy resolution of any existing 
telescope, and can be neglected. 

The flux expected in the "high" pair-annihilation region of the LHT model 3 is plotted 
in Fig. 1. For comparison, the flux due to BhBh^ 27, computed in [11], is also shown. 
Throughout the parameter space, the Z7 final state provides a stronger monochromatic 
photon flux, with the ratio of the Z7 to 77 flux varying between about 1.5 and 100. (The 
largest values of this ratio are obtained for M ~ 300 GeV, where the 77 cross section 
is supressed due to an accidental cancellation.) Similar results are obtained in the "low" 
pair-annihilation and the coannihilation regions of the parameter space. 

Experimental searches for the anomalous high-energy gamma rays using Atmospheric 
Cerenkov telescopes (ACTs), such as HESS [18] and VERITAS [19], are currently under 

3 In all calculations of photon and neutrino fluxes, we assume M = Mt_ — f in the pair-annihilation 
regions and M T _ = f in the coannihilation region. The impact of varying these parameters on the flux 
predictions is very small. 
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way, and the space-based GLAST telescope [20] is expected to perform such a search start- 
ing in 2007. Numerical simulations indicate that dark matter density may have a sharp 
peak at the galactic center, in which case the flux of the anomalous gamma rays would be 
maximized for a line of sight towards the center of the Milky Way. The principal source of 
background for a search focusing on the galactic center region is the recently discovered pow- 
erful point-like gamma ray source [21], whose spectrum strongly suggests that its nature is 
astrophysical. Taking into account this background, Zaharijas and Hooper [22] estimate that 
the minimal WIMP-related monochromatic photon flux required for 5<r discovery is about 
10~ n cm _2 sec _1 for GLAST (which is sensitive to photon energies up to about 300 GeV) 
and 10~ 12 cm^sec" 1 for the ACTs (sensitive to £ 7 > 200 GeV). The Z7 line predicted in 
the LHT model would be observable at GLAST and the ACTs only if the fluxes are enhanced 
by a strong spike in the dark matter concentration around the galactic center: depending 
on the Bh mass, values of JAQ in the 10 2 . . . 10 4 range are required. Many (though not all) 
models of the galactic halo contain such spikes: for example, the profile of Moore et al. [23] 
predicts J w 10 5 for AQ = 10 -3 sr. 

However, the prospects of the future searches are further restricted by the constraints on 
the continuous component of the photon flux from the observations of the gamma rays from 
the galactic center in the energy range up to 30 GeV by EGRET [24, 22]. The continuous 
component of the flux in the LHT model was computed in Ref. [11]. The maximum value of 
the Z7 photon flux compatible with the EGRET constraint (independently of the galactic 
halo profile) is shown in Fig. 2. It is clear that the maximum flux is substantially below the 
sensitivity of the current and near-future telescopes throughout the parameter space. 

Given the presence of a point-like background source at the galactic center, other regions 
with possible dark matter overdensity were suggested as potential targets for a search for 
anomalous gamma rays. These include dwarf spheroidal companion galaxies to the Milky 
Way such as Saggitarius, Draco and Canis Major [25, 26], as well as the Large Magellanic 
Cloud and the M87 galaxy [27]. In addition, in models where the galaxy is built up from 
hierarchical dark matter clustering one should expect localized clumps of dark matter inside 
the Milky Way halo. The values of JAQ expected for these objects are model-dependent. 
For example, the dark matter profiles in dwarf spheroidals surveyed in Ref. [25] give JAQ ps 
10~ 3 — 1 for AQ = 10~ 3 , 10~ 5 . For dark matter clumps in the halo, Baltz et. al. [29] estimate 
that a typical clump would have JAQ ps 0.4 at AQ = 1.5 x 10~ 4 . (This estimate uses 
the data from a simulation by Taylor and Babul [28].) For a space-based telescope such 
as GLAST, the background flux for these targets can be estimated by a simple power-law 
extrapolation of the extragalactic gamma ray flux measured by EGRET [30]: 



where k = 8.2 x 10~ n cm~ 2 s _1 sr _1 GeV _1 . Assuming a telescope with energy resolution 
SE/E = 0.1, the signal/background ratio is close to 1 for sources with JAQ ~ 1 in a model 
with M ~ 100 GeV. (The signal/background ratio decreases with increasing M and/or 
decreasing J.) However, the number of Z7 events expected at GLAST (A = 10 4 cm 2 ) is well 




2.1 



(14) 



6 



-10 

-11 

-12 

Log 10 5> ( s" 1 cm -2 ) 

-13 
-14 

150 200 250 300 350 400 450 500 
M (GeV) 

Figure 2: Maximum photon fluxes consistent with the EGRET bound on the continuous 
photon flux from WIMP annihilation. Solid line: "high" pair-annihilation region; dashed 
line: coannihilation region with = 300 GeV. Also shown are the sensitivities of GLAST 
(dotted red line) and HESS (dotted blue line) [22]. 



below 1 event/year, so that no discovery is possible. Of course, this pessimistic prediction 
could be proven wrong if dark matter turned out to be significantly stronger clumped at 
short scales than presently thought, resulting in larger values of J. Barring this possibility, 
a telescope with a larger effective area (A ^ 10 7 cm 2 ) would be required to begin probing 
the heavy photon dark matter model in this channel. 

The final annihilation process giving rise to monochromatic photon flux, BhBh—> hj, is 
strongly suppressed. It cannot occur via an s-channel Higgs exchange, because the process 
h* — > hj is forbidden by the Ward identity of Quantum Electrodynamics. The proof is as 
follows: Let the photon and on-shell Higgs have momenta k and p, respectively. The matrix 
element has the form 

M (h* -> hi) = e* (k) M" , (15) 
where Xi^ can be decomposed as 

M* i = A (k, p) k^ + B (k, p) p M . (16) 

The Ward identity requires that k^M.^ 1 = 0. For the first term this is true since k 2 = 0. For 
the second term, p-k = only if the initial Higgs is on-shell (when the process is kinematically 
forbidden); thus, B = 0. However, by the transverse nature of the polarization vector the 
first term provides no contribution to the matrix element; therefore, M. (h* — > hj) =0. As a 
result, BhBh^ h^y may only proceed via box diagrams with T-odd and T-even fermions in 
the loop, which receive no resonant enhancement. Moreover, the BuQq and BhLI couplings 
are of order <//10 ~ 0.035, further suppressing the cross section. 
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4 Neutrino Fluxes from Bh Annihilations in the Sun 
and the Earth 



Neutrinos produced in annihilations of heavy photons collected in the gravitational wells of 
the Sun and the Earth provide another potentially observable indirect signature of Little 
Higgs dark matter. The procedure for evaluating the neutrino fluxes in a given model is well 
established; a thorough review (in the context of SUSY) is given in [31]. Here we will follow 
this procedure to compute the neutrino fluxes expected in the LHT model. 
The number of WIMPs N collected in the Sun or the Earth obeys 

N = C-AN 2 , (17) 

where C is the capture rate and A is the annihilation rate per WIMP. Schematically, the 
capture rate is given by 

C ~ C T^~ (18) 
M 2 GeV 

where M Ge v is the WIMP mass in units of GeV, the quantity (cpb) is essentially the weighted 
average of the elastic WIMP-nucleus scattering cross sections (in pb) over the atomic com- 
position of the Sun or the Earth, and c is a coefficient determined by the properties of the 
astronomical body in question: c ~ 10 30 s" 1 for the Sun and c ~ 10 20 s _1 for the Earth. (See 
Ref. [31] for a more detailed discussion.) The only input from particle physics required to 
compute C is the elastic scattering cross sections, which were computed in Ref. [11] for the 
LHT model. The annihilation rate per WIMP A is schematically given by 

A=<^>, (19) 

where <r an is the total WIMP annihilation cross section, the average is over the thermal 
distribution of the WIMPs captured in the Sun or the Earth, and V eS is the effective volume 
of the "WIMP-sphere" inside the astronomical body. (For details, see Ref. [31].) 

Having computed C and A, we solve Eq. (17) to obtain the total WIMP annihilation 
rate: 

1 



- C tanh 2 (t VAC) , (20) 

where t ~ 4.5 • 10 9 years is the age of the Solar System. The experimental technique best 
suited to searching for high-energy neutrinos from WIMP annihilation relies on observing 
an upward-going muon created by a charged- current interaction of a muon neutrino in the 
rock below the detector. The rate of such muons per unit detector area is given by [31] 

Tdctcct = (2.54xl0- 29 m- 2 yr- 1 )%M 2 cV ^aAE^( iV ^) P .( M ) ) ( 21 ) 

s i F 

where % are the possible neutrino types, Oj and 6j are (known) coefficients describing the 
neutrino scattering and muon propagation in the rock, and F are the possible final states 
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Figure 3: The rate of neutrino-induced upward-going muon events expected from the heavy 
photon annihilation in the Sun. Solid line: "low" pair-annihilation region (the flux in the 
"high" pair-annihilation region is similar); dashed line: coannihilation region with = 300 
GeV. Also shown is the expected sensitivity of the IceCube detector (red/dotted line). 



of WIMP annihilation with branching fractions Bp. In the LHT model, the dominant 
annihilation channels are W + W~ and ZZ. The quantity 



is the second moment of the spectrum of neutrino type i from final state F. Here E t h is 
the threshold energy of the detector, and dNp^/dz is the neutrino spectrum, normalized per 
single WIMP annihilation into the final state F. This spectrum is a convolution of the initial 
neutrino spectrum at the production point with the propagation effects (including neutrino 
oscillations and absorbtion) on the way to the detector. In this analysis we used the neutrino 
spectra computed by Cirelli et. al. [32] , 4 and assumed a detector with a threshold energy 
of 50 GeV, representative of the IceCube experiment [33]. Note however that lowering 
this threshold would not have a substantial effect on the rates, since the sub-threshold 
contribution to the rate scales as (E t h/M) 3 and is at most of order 10% throughout the 
interesting parameter range in the LHT model. 

The rate of neutrino-induced upward-going muon events expected from the heavy photon 
annihilation in the Sun is shown in Fig. 3. The expected sensitivity of the IceCube detector 
is shown for comparison. The maximum possible rate (achieved at the low end of the allowed 
LTP mass range, M « 100 GeV) is about 1 event/yr/km 2 in the coannihilation region and 
about 0.5 events/yr/km 2 in pair-annihilation regions. Unfortunately these rates are well 
below the sensitivity of the IceCube. The sensitivity would need to be improved by a factor 
of a few hundred to a thousand before the fluxes predicted in the LHT models can be probed. 

4 We thank Marco Cirelli for providing us with the updated version of the spectra. 
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The rates of events due to heavy photon annihilation in the Earth are even smaller, of the 
order 10~ 5 events/yr/km 2 or below throughout the parameter space. 

5 Conclusions 

In this paper, we discussed two signatures of the heavy photon dark matter predicted by 
the LHT theory. We computed the cross section of the process B H B H — > Z7, leading to a 
monochromatic gamma ray signature. We also computed the flux of high-energy neutrinos 
from the annihilation of the heavy photons trapped in the Sun and the Earth. Unfortunately, 
the near-term prospects for observing both signatures are rather poor, the gamma ray 
signature could in principle be observed by GLAST, but only if dark matter is very strongly 
clumped at short distance scales. In the neutrino case, the predicted flux is too small to be 
observed at the IceCube. 

This study complements two previous analyses of the discovery prospects for the LHT 
dark matter. While all predictions are subject to significant astrophysical uncertainties, it 
appears that the most promising search channels are the "secondary" gamma rays produced 
in hadronization and fragmentation of the primary WIMP annihilation products [11], and 
anomalous high-energy positrons [12]. In these channels, the signal may be observed by the 
near-future instruments, GLAST in the case of gamma rays and PAMELA and AMS-02 
in the case of positrons. On the other hand, the LHT model predicts that no signal will 
be observed in the near-future direct detection and high-energy neutrino searches, while 
the monochromatic gamma ray signal is very unlikely. By testing the predicted pattern of 
signals, astroparticle experiments will provide an important test of the LHT dark matter 
hypothesis, complementary to the more direct searches for the new particles predicted by 
the LHT model at the Large Hadron Collider. 

Acknowledgments — We are grateful to Dan Hooper for helpful correspondence, and to 
Andrew Noble for useful discussions. We thank Marco Cirelli for clarifying the results of 
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